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Abstract 

A new formula for the gravitational potential of flattened systems is proposed. It is a modification of the 
Miyamoto-Nagai potential and should be applied to very flattened systems, exponential discs as a typical 
example. The resulting rotation curve agrees sufficiently well with that obtained by using special functions 
and the total masses remain the same. The functions contained in the new term can improve the agreement 
for the rotation curve and also reduce the effect of negative density values which appear off the midplane. 
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1 INTRODUCTION 

The gravitational potential proposed by Miyamoto and 
Nagai (1975) (MN) is well known. Its applications to 
models of the Milky Way and other disc galaxies are 
numerous (e. g. Allen Santillan 1991, Ninkovic 1992, 
Dinescu et al. 1999, MacMillan & Binney 2011). The 
observational evidence is in favour of exponential discs 
characterised by the dependence of the circular speed on 
the distance to the centre (rotation curve) as given in 
Freeman’s paper (1970). As to the fitting of the circular- 
speed dependence resulting from the Miyamoto-Nagai 
formula to Freeman’s curve, there are some difficul¬ 
ties indicated in an earlier paper of the present author 
(Ninkovic 2003a). In the present paper a new formula 
for the gravitational potential is proposed. It appears as 
a generalisation of MN and it contains the formula anal¬ 
ysed in the earlier paper (Ninkovic 2003a) as a special 
case. The density generating this potential as function 
of the coordinates is also studied. 


2 THE POTENTIAL FORMULA 

The gravitational potential in stellar systems (subsys¬ 
tems), such as, for instance bulges and discs of spiral 
galaxies, is often represented by a variant of the Green 
formula (e. g. Cuddeford 1993). This would mean that 
in general the gravitational potential $ should have the 
following form 


$ = 


GM 

~1T 


( 1 ) 


where G and M are constants, whereas TZ is a func¬ 
tion. The first constant (G) is the universal gravita¬ 
tion constant, the second one (M) is the total mass of 
the system (subsystem) appearing as the source of the 
gravitational field. 1Z is a function of the generalised 
coordinates and time. It must satisfy the following con¬ 
ditions: to have dimension of length and to be approx¬ 
imately equal to the distance from the centre of the 
source system r at the points lying very far from the 
system centre. 

The potential of a flattened stellar system, as here 
of interest, is assumed to depend on two arguments: R, 
distance to the axis of symmetry z, and \z\, distance to 
the midplane. The function proposed by Miyamoto and 
Nagai (1975) is a well-known example 


TZmn = 


R 2 + (a + v4 2 + 6 2 ) 


1/2 


( 2 ) 


In this function a and b are constants and for signifi¬ 
cantly flattened systems satisfy b/a < 1. When the func¬ 
tion TZmn is substituted in equation (1), one obtains the 
Miyamoto-Nagai potential formula. The present author 
proposes here a potential formula which contains a new 
term 1Z jv, as follows 


$ = 


GM 

TZmn ~ TZn 


( 3 ) 


In a general case the new term TZn is also a function 
of the same two variables, R and \z\ and is everywhere 
positive. Since TZmn satisfies the condition of being ap¬ 
proximately equal to the distance r at the points lying 
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very far from the system centre, the new term at such 
points must be negligibly small compared to TZmn■ In 
this way the difference TZmn — TZn will satisfy the con¬ 
dition of being approximately equal to the distance from 
the centre r at very distant points. 

3 THE ROTATION CURVE 

The surface density of the disc of a spiral galaxy, bear¬ 
ing in mind the observational constraints, is usually as¬ 
sumed to be exponential (often referred to as Freeman 
law), for instance its luminosity or mass density p obeys 
the following formula (e. g. Deg & Widrow 2013, then- 
equation (13)) 


is about 1.5 times as large as that of the exponential 
disc (e. g. Ninkovic 1992). 

The purpose of introducing the new term IZn (equa¬ 
tion (3)) is to reproduce the property of intersection 
with the Keplerian curve and to have the same to¬ 
tal mass as for the exponential disc. As shown earlier 
(Ninkovic 2003a), this is possible even with a constant 
substituted for TZn- The constant introduced here will 
be Rd, the scale length. Then in the fitting the rotation 
curve one should determine the ratio (a + b)/Rd- Ac¬ 
cording to the definition assumed in the present paper 
(equation ( 1 )) the potential cannot have negative val¬ 
ues, so this ratio must be greater than 1. We find a best 
fit value of 2.1. The fit can be further improved. This is 
done by generalising the term TZn in the following way 


p(R,z) = p(0)exp(^ - sech2 (^i^j > ( 4 ) 

where Rd and Zd are constants. As easily seen, the sur¬ 
face density following from this expression will depend 
on R as a simple exponential function with Rd as the 
scale length. 

Unfortunately, equation (4) yields no analytical solu¬ 
tions for the potential via the Poisson equation. This is 
important because with the known potential it is pos¬ 
sible to obtain the circular speed u c by using the well- 
known relation 


u c = 



z = 0 


(5) 


The plot of circular speed versus distance R according 
to Freeman’s solution (Fig. 1) is characterised by exis¬ 
tence of radius at which the circular speed is the same 
for the Keplerian and exponential disc cases of the same 
mass. As a consequence, at infinity the circular speed 
for the exponential disc approaches the Keplerian curve 
from ’’above” , i. e. from higher values of u c - This is 
not the case with the potential formula of Miyamoto 
and Nagai. Substituting equation (2) in equation (1) 
and applying (5) one obtains the dependence of the 
circular speed on R corresponding to the Miyamoto- 
Nagai potential. Since z = 0, the parameters a and b en¬ 
ter the circular-speed dependence always through their 
sum a + b. In other words, the flattening expressed by 
means of b/a, does not affect the circular-speed depen¬ 
dence. The maximum circular speed occurs at about 
1.4(a + 6 ) and the circular-speed value resulting from 
the Miyamoto-Nagai formula is everywhere smaller than 
the value yielded by the Keplerian dependence for the 
same total mass. Due to this a good fit between the cir¬ 
cular speed following from the potential of Miyamoto 
and Nagai and that corresponding to the exponential 
disc is achieved at the cost of different total masses. 
Usually the total mass in the Miyamoto-Nagai formula 


TZn 





( 6 ) 


where 71 < 0.5, 72 < 0.5 Though it may seem that the 
number of parameters now tends to be too large, this 
is not the case in practice. For instance, Rd is quite ac¬ 
ceptable to be substituted for ci. The second term in 
equation ( 6 ) (the function of |z|) is foreseen because of 
the density distribution the potential returns (see next 
section). With regard to all equations written above 
except (4) one obtains the following expression which 
yields the circular speed 


u c = 


GMR 


\l 


\/ R 2 + (« + b ) 2 - \Rd ( l +( l + T§) 71 ) 


R 


\ \JR 2 + (a + b ) 2 n Rd 


r r r 2 

71 D I + R 2 


71-1 


(7) 


The curves presented in Fig. 1 follow from equa¬ 
tion (7) with a + b = 2.1 Rd, 71 = —0.1 and 71 = —0.3. 
They are presented together with the corresponding 
curve from Freeman’s (1970) paper. The dimension¬ 
less variables, £ and y, are defined as: £ = R/Rd , X = 
y/GM/Rd- The comparison requires the scale lengths 
Rd and the total masses Ai to be equal in both for¬ 
malisms. 


4 THE DENSITY 

The density which generates a gravitational potential 
and the potential are related through the Poisson equa¬ 
tion 


V 2 $ = -4t xGp . ( 8 ) 
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Figure 1. Rotation curves from equation (7) as described in the 
text (blue curve 71 = —0.1 , red curve 71 = —0.3) together with 

Freeman’s (1970) curve; £ = R/Rd , X — u c 



Figure 2. Density dependence on height \z\, for R = 3 Rd, 71 = 
72 = 0 (see equation 6 ); distance unit Rd, density unit A4 R 


V 2 is the Laplacian, p is the density. The full formula 
is given in the Appendix. 

In the density calculation the constants Ci and C2 
(equation ( 6 )), as well as the ratio b/a, must be speci¬ 
fied. The fit for the rotation curve has yielded a + b = 
2.1R d , under the condition 6 Co (flattened mass distri¬ 
bution) , so it is easy to conclude that a is approximately 
2R d . In such conditions the first term within the brack¬ 
ets in equation ( 6 ) should not differ substantially from 
1 . As for the second one, a too small C2 can contribute 
to significant values of 7 Zn, for instance C2 = b, leads to 
such a high 7 Zn that the denominator in equation (3) 
becomes negative and, as a consequence, the potential, 
contrary to the convention assumed here, becomes also 
negative. Thus, it is reasonable to assume c± = C2 = R d . 

In spite of introducing the new term TZn (equations 
(3) and ( 6 )) negative density values cannot be avoided. 
In the midplane the density is nowhere negative. Its de¬ 
pendence on R in z = 0 is a monotonously decreasing 
function. As for the dependence p(z), R = const, there 
are three types of profiles: i) with a minimum p m i n < 0 
and then approaching zero from the negative side (Fig. 
2 ); ii) a wavy profile, with one minimum and one maxi¬ 
mum and then approaching zero, along R = 0 both ex¬ 
trema can be positive, in general the minimum is nega¬ 
tive, the maximum positive (Fig. 3); iii) a sound profile, 
without negative values and monotonously decreasing, 
but obtained along R = 0 only (Fig. 4). In order to be 
more clearly visible the most essential part of the curve 
is magnified and given in the same figure (Fig. 2 and 
Fig. 3). 

In the simplest case - 7i = 72 = 0 (equation ( 6 )) - the 
density dependence p(z), R = const, is always of type 
i) (Fig. 2). The z value at which the density minimum 
occurs depends on R, the smaller R is, the closer is the 


minimum to the midplane. Finally, when R tends to 
infinity, the z value also tends to infinity, because the 
density tends to zero. 

When the term 7 Zn (equations (3) and ( 6 )) takes a 
functional form, i. e. both 71 and 72 are non-zero, the 
dependence p(z), R = const, changes the form, towards 
the types ii) and iii). For 7 ! some kind of critical value 
appears at which the dependence p(z), R = 0, obtains 
the form iii) (Fig. 4). It should be 71 < 0, the modulus 
is affected by the ratio b/a. For instance, for the pair of 
values b = 0.18, a = 1.92, the critical value is 71 = —0.3, 
for b = 0.08, a = 2.02, the corresponding is 71 = —0.5, 
for b = 0.02, a = 2.08, 71 = —0.6, etc. The influence of 
72 is not so strong; it is better to be 72 > 0 , but increas¬ 
ing 72 leads to weaker decreasing of p with increasing z. 
The most important is that increasing the modulus of 

71 further does not improve the situation, as the nega¬ 
tive density values remain and the wavy profiles survive 
(Fig. 3). However, the moduli of both exponents 71 and 

72 are limited. As for the former, the main limitation 
comes from the fit of the rotation curve, in the case 
of the latter, since it is supposed to be positive, the 
physical limitation (72 < 0.5) becomes essential. The fi¬ 
nal conclusion is that the best achievable result is to 
obtain a monotonously decreasing density along R = 0 
and wavy density profiles along cylinders R > 0. For the 
purpose of giving a more clear explanation the follow¬ 
ing example is chosen. The values of the parameters in 
equations (3) and ( 6 ) are: a + b = 2AR d (b = 0.18R d ), 
ci = C 2 = Rd, 71 = —0.3, 72 = 0.3. The rotation curve 
(Fig. 1 - red curve) and the density profiles (Fig. 3 and 
Fig. 4) correspond to this example. 

In Fig. 2 we show another density comparison, typ¬ 
ical for 71 = 72 = 0 , corresponding to the same a + b 
and the ratio b/a. The wave (Fig. 3) can be explained 
by the rather rapid density decrease with increasing dis- 
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Figure 3. Density dependence on height \z\, for R = Rd, 71 = 
—0.3 ,72 = 0.3 (see equation 6 ); distance unit Rd , density unit 

MR,f 
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5 DISCUSSION AND CONCLUSIONS 

The present author proposes an analytical form for 
the gravitational potential which would correspond to 
the exponential disc, the most luminous subsystem of 
a spiral galaxy. This is done by modifying the well- 
known formula of Miyamoto and Nagai (equations (1)- 
(3) and ( 6 )). In this modification there are three essen¬ 
tial parameters: the total mass Ai, the exponential scale 
length Rd and the other scale length b. The other three 
parameters - a , 71 and 72 are auxiliary. Their role is to 
improve the fit of the rotation curve particularly, their 
values are limited. Since the system (subsystem) under 
study is flattened, it must satisfy a^$> b, and bearing in 
mind the constraint based on the rotation curve, one ob¬ 
tains a « 2 Rd- The exponents 71 and 72 cannot exceed 
0.5 and 71 is constrained to be small and negative, while 
72 lies in 0 < 72 < 0.5. Thus an exponential disc mod¬ 
eled in the way as proposed here should be generally 
characterised by its total mass and two scale lengths. 

Though negative density values cannot be avoided 
completely, the functions in the new term (equation ( 6 )) 
contribute to a significant reducing of this effect. A dis¬ 
cussion concerning the analogous spherical-symmetry 
model can be found in Ninkovic, 2003b. 
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A APPENDIX 

Equations (1) and (3) from the main text are given again 


Figure 4. Density dependence on height \z\, for R = 0, 71 = —0.3 
, 72 = 0.3 (see equation 6 ); distance unit Rd, density unit A4 R^ 3 



tance to the midplane. Beyond the centre the density 
values in the midplane are generally low, significantly 
lower than at the centre, so that at reasonable distances 
to the midplane they are practically zero. Here one ob¬ 
tains the density values by using the Poisson equation 
(8) which in the case of axial symmetry means as the 
algebraic sum of three terms (see Appendix). Clearly, 
in a numerical procedure instead of density values of 
exactly zero we usually have values of low moduli with 
both signs. In this way the wavy density profiles ob¬ 
tained here can be easily understood as due to numeri¬ 
cal precision issues. 
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TZ = IZmn — TZn ■ 

The designations used are: 'll, the gravitational potential, 
G the universal gravitation constant, whereas IZmn and IZi v 
are two functions depending on the arguments R and 2 (axial 
symmetry). 

The Poisson equation ((8) from the main text) relating 
the density and potentail is also rewritten 

V 2 $ = -4t rGp . 

With regard to the axial symmetry the Laplacian has the 
form 

<9 2 $ 1 d<5? d 2 <f> 

V $ =-1-1-. 

dR 2 RdR dz 2 
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The derivatives will be: 


<9$ 

dR 


GM 
IZ 2 


dTU 


dR 


dIZ N 

OR 


81Zn z ( z 

~dT^ l 2 lc\ + r ! 


2 \ 72-1 


d 2 $ 

dR 2 


2 GM 
TO 


GM 


&1Zmn DIZn 

dR 8R~ 


8 2 TZmn 8 2 1Zn 


iz 2 


d 2 <P _ 2 GM 

dz 2 ~ R? 

GM 

~ iz 2 


8R 2 


8R 2 


cHZmn cHZn 


dz 


dz 


d-llMN 9 “IZn 
dz 2 dz 2 


The first function in the denominator IZmn and its nec¬ 
essary derivatives are: 


IZmn = 


R 2 + (a + Vz 2 +b 2 Y 


dIZ M n R 


1/2 


dR IZmn 


8 2 IZmn _ Rmn ti M j, 


d 2 n N n , ^ z 2 A z 2 \ 72-2 

~d^~ - 2 ^ 2 ( 72 - 1)^3 { 1+ WJ + 

R*\ Rl, 
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The second function in the denominator IZn and its nec¬ 
essary derivatives are: 


IZn = ~Rd 


, -R V 1 (, * 

Rl) + \ Rl 


dIZ N R /, J? 2X71_1 
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